function parn=clean_gould_par(par)
% parn=clean_gould_par(par) computes mean bias,
% adds bias to par series and sets parn<0 to
% 0 to simplify plotting. 
% 4/10/01

m=find(par<=0);
parb=mean(par(m));
parn=par-parb;
m=find(parn<0);
parn(m)=0*parn(m);